The effects of stochastic population dynamics on food web structure 
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We develop a stochastic, individual-based model for food web simulation which in the large- 
population limit reduces to the well-studied Webworld model, which has been used to successfully 
construct model food webs with several realistic features. We demonstrate that an almost exact 
match is found between the population dynamics in fixed food webs, and that the demographic 
fluctuations have systematic effects when the new model is used to construct food webs due to the 
presence of species with small populations. 
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I. INTRODUCTION 



There is a crucial distinction in biology between the 
proximal causes of phenomena and the evolutionary ex- 
planation [l| . While the former is useful in giving specific 
insight into a system being studied, it is only through the 
understanding of the evolutionary mechanisms in action 
that general principles can be understood. Applying this 
understanding to ecosystems, the proximal explanation 
of food web structure relates to the identity and habits 
of the species within it [for example 0, H[ , but the un- 
derlying, evolutionary explanation requires understand- 
ing of the emergence and stability of those inter-specific 
relations which contribute to the food web. Whereas al- 
ternative models of food web structure such as the niche 
model [H and the cascade model [f| are essentially static 
with respect to their species composition, the Webworld 
model [6[ has been used to demonstrate the formation of 
complex food webs stable to evolutionary perturbation 
with several emergent features consistent with observed 
ecosystems @, BUI- O ne aspect of real ecosystems which 
it has not been possible to reproduce in the Webworld 
model is the existence of oscillations about a fixed point 
in the population dynamics. Historic studies of popula- 
tion dynamics have tended to focus on systems which ex- 
perience oscillations, such as the Snowshoe Hare [l(|, or 
Great Tit ll| , but the Webworld model shows monotonic 
approach towards a fixed point. In part this is a conse- 
quence of using population dynamics continuous in time, 
since chaotic dynamics often emerge when large over- 
shoots of the carrying capacity are possible [l2| , an effect 
mitigated by using very small time steps or effectively 
continuous dynamics. Kfivan & Eisner [131 ] demonstrate 
the occurrence of limit cycles and chaotic behaviour even 
using continuous population dynamics and a Holling type 
II functional response. It is likely that the explicit compe- 
tition term in the functional response used by Webworld 
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prevents these cycles occurring. It was, however, demon- 
strated by McKane & Newman [l4[ that cyclic behaviour 
could arise in a model using finite populations which dis- 
appears when the infinite population limit is taken. Since 
the continuous population model of Webworld is equiv- 
alent to this limit, one explanation for the absence of 
realistic population cycles is the absence of individuals 
from the model, whose demographic stochasticity could 
potentially lead to resonant oscillations in the population 
dynamics. 

Stochasticity in ecosystem models tends to be divided 
into the three categories described by Lande et al. [H|]. 
The first of these, demographic stochasticity, is simply 
the result of considering small numbers of entities in a 
system. For large populations, the law of large num- 
bers tends to ensure that birth and death processes, 
happening at random to unrelated individuals, never- 
theless average out so that differential equations are an 
adequate description of the processes. For small pop- 
ulations it may be essential to the future of a species 
how many, when, and in what order births and deaths 
occur. In the case of most ecosystems, including those 
modelled by Webworld, there exist at least some species 
with small populations If these species exert 

an influence on the ecosystem as a whole, for instance 
by having a role as a top predator [171 ] , the dynamics of 
the system may crucially depend on the demographics 
of species with small populations, and as a consequence 
demographic stochasticity cannot be neglected even for 
the understanding of the broader ecosystem. The re- 
lated effect of demographic heterogeneity, characterised 
by Melbourne & Hastings [l|| , concerns the demographic 
stochasticity within species of groups of heterogeneous 
individuals, for example males and females. Even if the 
population of the species as a whole were constant, the 
sex ratio is generally not strongly constrained by effects 
such as prey abundance in the same manner as the pop- 
ulation of the species as a whole, and hence is more vul- 
nerable to large-amplitude fluctuations. The Webworld 
model does not explicitly incorporate sexual reproduc- 
tion, or distinguish between individuals by sex, but in 
Section IIII CI we show that demographic heterogeneity 
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in terms of foraging strategy can have important con- 
sequences for the food web. The two further types of 
stochasticity discussed by Lande et al. [HJ are not exam- 
ined in this paper. Of these, environmental stochasticity 
and indeed environmental periodicity, could be incorpo- 
rated into the Webworld model, but the work of McKane 
& Newman indicates that the intrinsic fluctuations 
of demographic stochasticity are sufficient to have large 
scale consequences on the population dynamics. Because 
demographic stochasticity can be naturally generated by 
the probabilistic occurrence of events with characteris- 
tic frequencies, we feel that it is more natural to intro- 
duce this as the first form of stochastic influence on food 
web structure than to apply an external forcing, which 
would require additional assumptions about the level at 
which it acts (basal resource availability, death rate, etc.) 
and the degree to which the effects on different species 
were correlated. In this paper we therefore develop an 
individual-based model designed to replicate Webworld 
in the large-population limit. The final common type 
of stochasticity relates to measurement error. As a mea- 
surement of a real ecosystem will not yield an exact count 
of all individuals of each species, the dynamics of the real 
system will differ from those that would occur if the mea- 
surements were correct. As such, even measurements of 
a constant population will show apparent fluctuations. 
Since the Webworld mathematical model utilises exact 
populations, the population records are at any moment 
in time precise, and this form of stochasticity is not rel- 
evant in this case. 

A second aim of this work is to introduce a model to 
which the results previously published for Webworld are 
an approximation, valid under conditions which need to 
be properly understood. It is more difficult to obtain 
precise results for the more detailed model developed in 
this paper than for the parent Webworld model, either 
analytically or numerically. The parent model remains 
the most appropriate level of detail at which to derive 
gross properties of large ecosystems with relative compu- 
tational economy. In addition, food webs established in 
the parent model can be used to grow food webs adapted 
to stochastic population dynamics. Imported food webs 
can either be used as the basis for on-going evolution, 
or species can be drawn from a parent-model 'mainland' 
to construct immigrant communities, as is done in the 
present paper. By working towards a model in which 
individual actions can be described, the consequences 
of differing plausible behaviours can be examined. In 
particular, an individual-based model such as that intro- 
duced in this paper is the appropriate level of detail at 
which to examine foraging strategy, or can be considered 
as an essential step in deriving an agent-based model in 
which experience-derived strategy optimisation, and life- 
history effects, can be examined. The approach taken 
in this paper is to model species by discrete individuals, 
who are affected by events such as death or the slow di- 
gestion of consumed prey, and who interact in a pair- wise 
manner between predators and prey. This approach has 



the virtue that it can be formally described in the same 
manner as a set of chemical reactions, and the Gillespie 
algorithm 19] can be used to simulate essentially exact 
and typical histories. In Section [TT] we expand on the na- 
ture of this description; in Section III Al we recap those 
aspects of the Webworld model, more fully expounded 
in Drossel et al. [6|], which are necessary to derive the 
individual-based model. This derivation is presented in 
Sections III Bl to Hi Dl In Section Hi El we describe variants 
on the basic model relating to foraging strategy selection 
in order to investigate the sensitivity of the Webworld 
model to assumptions about decision making which can- 
not be adequately investigated using continuous popula- 
tion dynamics. A summary of the stochastic model is 
given in Section Hi Fl 

II. MODEL 

The Webworld model from which results have been 
published in previous papers @, H, 0, US, Hl| is stochastic 
to the extent that species are added as random mutants 
of a randomly selected parent species, but the behaviour 
of the system between speciation events is entirely deter- 
ministic. In particular, the population dynamics is repre- 
sented by a set of differential equations linking the rate of 
change of the population of each species to its own popu- 
lation as well as the populations of its prey, competitors 
and predators. This is notably similar to the represen- 
tation of the dynamics of a chemical system in terms of 
reaction rates. Using the insight that such deterministic 
reaction rate equations emerge from large numbers of in- 
teractions of discrete molecules, and the desire to model 
ecosystem dynamics in terms of individual organisms, in 
this paper we deduce the corresponding set of reactions 
which gives rise to the Webworld model as its expected 
behaviour. Interactions changing the state of individuals 
will be written after the manner of chemical reactions, 
which follow the general model 

Xa+Xb k ^ X c , [I] 

where X A and Xg are individuals present before the in- 
teraction, of types A and B respectively, and Xc is the 
individual of type C which remains after the interaction. 
Writing the abundance of these types in turn as Na and 
Nb, the total rate at which [I] occurs is given by 

R = kN A N B , (1) 

where k, written over the reaction arrow of [I], is the reac- 
tion rate coefficient. Concrete examples of reactions are 
given in Section III Bl In the Webworld model the rate of 
predation is determined by a ratio-dependent functional 
response, which we expect to be encoded in the appro- 
priate reaction rate coefficients. 

Having chosen to make the stochastic Webworld model 
reaction-based, an efficient means of implementing the 
dynamics in an essentially exact manner is to use the 
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Gillespie algorithm [19]. One of the assumptions we re- 
tain in order to use this method is that spatial effects are 
not relevant; in a chemical reaction scheme this is equiv- 
alent to supposing the system to be well-mixed. Under 
this assumption the reaction rates can at all times be cal- 
culated in terms of the abundance of the reactants with- 
out regard to their distribution. The original Webworld 
model makes similar assumptions about the absence of 
spatial structure and hence is well-matched. A model in- 
corporating spatial structure is likely to use spatial con- 
siderations as the basis of the functional response under- 
lying Webworld, and may therefore need to describe in 
detail the individuals of each species and their particular 
behaviour. In addition to the direct interest in the effects 
of demographic stochasticity on the assembly and func- 
tioning of ecological communities, the model described in 
this paper can be understood as a necessary step towards 
an agent-based model which could explore the effects of 
individual-level behaviour and spatial inhomogencity. 

For conciseness we will refer to the Webworld model 
with deterministic population dynamics, discussed in 
Drossel et al. [6( and briefly introduced in Section Hi Al as 
the parent Webworld model, and the model variant de- 
veloped in this paper as the stochastic Webworld model. 
The term Webworld by itself refers to features generic to 
both. 



A. Webworld 

The Webworld model consists of dynamics on three 
separated time-scales, of which the longest corresponds 
to the introduction of new species. A detailed descrip- 
tion of the parent model can be found in Drossel et al. 
6] . In the present paper we do not consider the introduc- 
tion of species through evolution, but instead examine 
two types of system. In the first, an existing commu- 
nity, constructed in the parent model, is interpreted in 
terms of discrete populations, and the evolution of the 
stochastic population dynamics is examined. In the sec- 
ond type of numerical experiment, we use a large com- 
munity assembled in the parent model as a species pool, 
from which individuals are drawn during the construction 
of an immigrant community, in which the population dy- 
namics are again stochastic. The subsequent discussion 
of the parent model will only restate those features of 
the parent model which are necessary to construct this 
non-evolutionary stochastic model. 

Interactions between species in the Webworld model 
are governed by two evolved matrices, scores S and com- 
petition a. Each pair of species i,j has an associated 
score Sij — —Sji describing the ability of one species 
to feed on the other. A positive value of Sij indicates 
that species i is able to feed on species j, larger values 
corresponding to preferred prey. The anti-symmetric na- 
ture of the matrix S thus precludes mutual feeding and 
cannibalism. The competition matrix, a, describes the 
degree to which species compete with one another when 



feeding on the same food sources; this is described by 
a symmetric matrix, which has a minimum value of 0.5 
for highly distinct species and a maximum value of 1 
for infra-specific competition. In the parent Webworld 
model S and a are derived by an intricate system of evo- 
lution, but it has also been demonstrated [221 ] that varied 
communities can be assembled by immigration from an 
evolved species pool. 

The intermediate timescale of the Webworld model cor- 
responds to the population dynamics. By encapsulating 
the per-capita income of species i from prey species j in 
the functional response, gij, the rate at which the popu- 
lation of each species i changes is described by 



N, 



= A ^ 9ijNi - 9jiNj - '/A , 



(2) 



where the first term corresponds to the gain of resources 
through feeding, the second term corresponds to losses to 
predation, and the final term relates to the death of indi- 
viduals through other processes. Webworld equates en- 
ergy resources to population, equivalent to an assumption 
of equal body size for all individuals. In these terms, the 
ecological efficiency, A = 0.1, which appears in the first 
term of ([2]), corresponds to the fraction of prey resources 
which can be converted into new predator individuals. 

The shortest time-scale of the Webworld model relates 
to the adaptation of the functional response to changes 
in population. Part of this response is a direct conse- 
quence of the reduction in the prey density, which re- 
duces the rate at which individuals can be found. The 
functional response also incorporates adaptation by for- 
agers, who select prey to maximize their income. It is as- 
sumed that species as a whole rapidly respond to changes 
in the relative abundance of their prey and of their com- 
petitors. Section III El discusses an alternative method 
we consider by which the evolutionarily stable foraging 
strategy can be identified, but an ecologically plausible 
mechanism is likely to require an agent-based and spatial 
model, which will provide a mechanism capable of repro- 
ducing the functional response. The exact form of the 
functional response as used by Webworld is 



9ij 



fij Sij Nj 



bN j +J2k ^ikfkjSkjNk 



(3) 



which uses the scores, S, and competition, a, already de- 
scribed. The term b in the denominator limits population 
growth for predator i of species j when Nj 3> Ni, whereas 
the sum limits the population of competing predators 
when their population is non-negligible. The term fij 
represents the fraction of its time that predator i spends 
feeding on prey species j, and was demonstrated by 
Drossel et al. [y] to have an evolutionarily stable strategy 
given by 



fi. 



9ij 



Efe 9ik ' 



(4) 



which can be found by iteratively solving ^ and 
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B. Target reaction scheme 

An essential modification made to the population dy- 
namics in creating this stochastic model is the division 
of each species into hungry individuals, who are actively 
involved in feeding, and sated individuals, who are not. 
This standard assumption is used to derive a functional 
response in which the predation rate for high prey abun- 
dance has a limiting value [for example[23j . An important 
consequence of this division is that the rate of predation 
is not dependent on the total abundance of the predator, 
but of its hungry members. As prey become more numer- 
ous, it should be expected that a smaller fraction of the 
predator individuals remain hungry and the predation 
rate per predator individual becomes restricted. This 
behaviour replicates an essential feature of the Webworld 
functional response. 

Having identified the need to divide species into hungry 
and sated individuals, it follows from ([2]) that four generic 
reactions are required. The first term in ([2]) corresponds 
to feeding, and hence to the reaction 



x; + x r 



[ii] 



where X[ and X* indicate an hungry and a sated individ- 
ual of species i respectively, Xj is any member of species 
j, and reaction rate coefficient k\ remains to be deduced. 
To restore individuals to the hungry state it is necessary 
to introduce the reaction 



X* 



[III] 



since no interaction is expected to be necessary for sated 
individuals to become hungry. To enforce the principle 
of energy conservation not only on average, but in the 
details of the reaction scheme, reproduction can be sup- 
posed to occur only after feeding, and to be limited to one 
reproduction event per prey individual consumed. This 
suggests an analogy to [In] , 



X* 



^2X 



[IV] 



where a sated individual of species i produces a hungry 
individual while becoming hungry itself. To complete 
the set of processes associated with ([2]) reactions corre- 
sponding to death are required. If death is equally likely 
to occur whether an individual is hungry or sated, the 
reaction can be written as 



Xi- 



[V] 



noting that, as with the prey individual in [IT], [V] applies 
to both hungry and sated individuals. 



of fixed abundance R. The expected rate of change of 
abundance of hungry individuals can then be written 



dN,, 



(5) 



where the four terms correspond to the four reactions 
in Section III Bl in the order they appear. iVj< and Ni* 
are the hungry and sated populations respectively, which 
sum to give the total population, iVj. Correspondingly, 
the number of sated individuals changes as 



N* = k x N %l R ~ k 2 N t , - h 3 Ni* - dNi. 



(6) 



with the population of the species as a whole changing 
according to the sum of ([5]) and ©, 



Ni = k s Ni* - dNi 



(7) 



The simplified version of @ appropriate to this scenario, 
using for the functional response, is 



Ni = XSl ° R N ~ dN 
bR + SioN 



(8) 



where the fixed resources have been identified as species 
0, and / = 1 occurs since no choice of prey is available. 
As described in Section fll A[ a = 1 for all intra-specific 
competition. 

Progress can be made by finding an appropriate steady 
state by which to eliminate time derivatives. By assum- 
ing the population of full individuals, N* , to be a con- 
stant, ([6]) becomes, 



Ni. = 



k x R + k 2 +k 3 + d 



(9) 



allowing us to eliminate from ([7]) the population of sated 
individuals. This leaves 



N 



k x k z R 



hR + k 2 



-N - dN. 



(10) 



which has obvious similarities to ([8]). Assuming that the 
k values do not depend on R, it follows that 



Losses to the resources through [II] occur at rate 



R 



(11) 



(12) 



and according to the Webworld model at rate S^R/ (bR+ 
SioN). In the limit of large R it follows that 



k- ? . 



S i0 - bd 



(13) 



C. Minimal reaction scheme 

To help deduce the reaction rate coefficients for the 
reactions in Section III Bl this section considers the sim- 
ple scenario of a single species, i, feeding on a resource 



Because SV,- is typically of order 10 for active feeding 
links whilst bd = 1/200, we neglect the latter term in the 
numerator of (1131) to obtain 



(1-A) 



Sio 



(14) 
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The fact that the rate coefficients of reactions . |IV| 
and [V] do not depend on any population implies that 
sated individuals can be considered as non-interacting 
entities. An implication of this in terms of prcdation is 
discussed in Section III Fl Of all occurrences of reactions 
and |IVj . a fraction A result in reproduction, in line 
with the ecological efficiency of the parent model. The 
reaction rate of the predation reaction is not so trivial, 
being given by 



h = 



1 

W 



(15) 



intra-specific competition causes the feeding rate of each 
member of the species to decrease in proportion to the 
total population, not merely the hungry population who 
are concurrently foraging. Section l!! Dl examines how this 
rate generalises to include inter-specific competition, and 
the effect of a diversified foraging strategy in reducing the 
intra-specific competition implied by (TH 



D. Foraging strategy and competition 

The generalisation of the reaction rate of predation 
must consider the effects of intra- and inter-specific com- 
petition. In the Webworld model, species adopting a 
mixed foraging strategy increase their reproduction rates 
by utilising more food sources, which is mathematically 
described by the reduction in intra-specific competition. 
Multiple predator species feeding on one prey are taken 
to be more effective at extracting resources from that 
prey species than a single predator species due to dif- 
ferences in strategy, and hence inter-specific competition 
is somewhat less detrimental than intra-specific competi- 
tion. The effect of a mixed strategy by a single predator 
species will be considered first. The optimal strategy, 
/, used by Webworld is the evolutionarily stable strat- 
egy (ESS) ; an individual cannot increase its reproductive 
success by adopting any other foraging strategy. Given 
the implication of the analysis in Section III CI that the 
probability of reproduction by a sated individual is A, 
and hence independent of its food source, it follows that 
the optimal strategy for a hungry individual also cor- 
responds to minimising the time before next becoming 
hungry. In general this time depends on the prey species 
chosen by an individual, and we write as the expected 
time before a hungry individual of species i, following a 
pure strategy of feeding on species j, next becomes hun- 
gry. For an optimal strategy ry is independent of j, and 
hence no individual can improve its reproduction rate by 
changing strategy. Importantly, this condition is not the 
same as minimising the time before next feeding, since 
the expected time in the sated state is 



(16) 



which is smaller for prey species corresponding to large 



Inspection of the Webworld functional response, ([3]), 
demonstrates that a single species following a mixed 
strategy consisting of two prey has instantaneous pop- 
ulation growth equivalent to considering two populations 
following pure strategies, where the sub-populations are 
given by = fijNi. In the absence of both inter- 
specific competition and predation, and denoting the 
prey as species 1 and 2, the total population growth is 
given by 



+ 



XS l2 N 2 
bN 2 + S i2 N i2 



N t2 = 
N l2 - 



\S iX Ni 



bNi + S a Na 
dNn-dN i2 , 



(17) 



which is equivalent to treating components Nn and N 2 
as separate species following pure strategies. In terms 
of the stochastic model, it follows that Nn is the total 
number of individuals either foraging for prey 1 or sated 
having eaten that prey, and hence we can write 



fijNi = fcjNv + N, 



(18) 



where <fi is the foraging strategy of hungry individuals 
and Ni*j is the number of individuals of species i sated 



by having consumed a member of species j. The depen- 
dence of (fT5|) on N*j can be eliminated by considering 
the balance of reactions affecting the population of sated 
individuals in the steady state, where 



fijNi 



,N, = ^-N, 



3 



dN* 



(19) 



The left hand side of (|19p represents the rate at which 
foraging individuals become sated, the terms on the r.h.s 
being the sum of hungering and reproduction, and death. 
Using (fl"8"|) to eliminate Ni*j from ([19]) . the foraging strat- 
egy of hungry individuals is found to be given by 



Ni_ 

Nil 



1 



bN, 



{Sij + bd) f^N 



f 



(20) 



We note that, subject to the condition that Sij ^> bd, 
(|2TJ]) can be written as 



fij Nj 



Nf 



bNj + fij Sij N Ni>N 



' fk 



(21) 



where the first fraction is exactly the form of gij in the 
parent Webworld model. Thus, in the absence of inter- 
specific competition, the rate of feeding through [IT] is 
given by 



<t>i<iNi>Ni 
hNvNj = 9t > ' J = gijNi. 

Jij^i 



(22) 



The appropriate foraging strategy for hungry individuals 
can be calculated using ([20]) . and (|22[) demonstrates that 
this results in dynamics consistent with the parent model. 

The reaction rate given by (|2"2"|) needs to be modified 
in the presence of inter-specific competition. By noting 
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the argument associated with (JTTj) , that species pursuing 
a mixed strategy can without loss of generality be split 
into species of pure strategy, it is possible to deduce the 
form of (f2"2"| for competition without strategy, and hence 
infer the behaviour of the system when both competition 
and strategy are allowed. First we equate the population 
growth equations of the parent and stochastic Webworld 
models for species i feeding on prey species j, 



XSijNi*j 



dN, 



XSijNjNiU 



bN j + J2k oiikfkjSkjNk 



dNj. 



(23) 

this can easily be simplified because the death terms and 
ecological efficiency appear equivalently on both sides. In 
steady state the number of sated individuals, Ni*j, can 
be found by setting the population growth of strategy i, j 
to zero, and hence 



kijifijNi-N^Nj 



■% + bd N 



■3' 



(24) 



where ku is the rate coefficient for the feeding reaction, 
[IT] , and the r.h.s. is the total rate at which the number of 
sated individuals is reduced. Rearranging in terms 
of the fraction Npj/ fijN gives 



becoming hungry, when reproduction might occur. For 
species i feeding on prey species j this time is 



J2k a ikfkjSkjNk b 



S l0 N 3 



Sij 



(30) 



where the first term is the reciprocal rate of feeding, and 
the second is the reciprocal of the total rate of hungering 
and reproduction. Given the definition of g in ([3]), it is 
clear that 



fij 

9ij' 



(31) 



For the ESS, the foraging strategy / is related to the 
functional response g by 



9-U 



J2k 9ik ' 



(32) 



and hence the expected time for an individual of species 
i to next become hungry is given by 



(33) 



kijNj = tk+} d fh N * 



- 1 



(25) 



which is independent of prey as expected. 



This ratio can also be obtained from (|23l) , from which we 
obtain 



fijNi _ bNj +J2k a ikfkjSkjN k 



N*j 



Substituting (26) into J25]) gives 



Sn + bd 



J2k a ikfkjSkjNk ' 



(26) 



(27) 



which is similar in form to ((3|). If we again assume that 
Sij 3> bd, this can be written 



^a ik f kj ^-Nk^ , 



(28) 



and is clearly seen to simplify to (| 1 5|) for a pure strat- 
egy in the absence of inter-specific competition, where 
an = 1. For the case in which a single predator has 
reduced intra-specific competition due to adoption of a 
mixed strategy, (f2"8")l gives 



1 



h ' J hjN 



(29) 



where the appropriate forager and prey populations are 
cf>ijNi' and Nj respectively, giving ([22]) . 

It is now possible to demonstrate that the ESS of the 
parent model corresponds as expected to the optimal in- 
dividual behaviour, which minimises the time until next 



E. Variants 

Two basic variations of the stochastic Webworld model 
are suggested by the details of the behaviour seen above, 
which relate to anti-predator behaviour and to decision 
making. Having shown that the rate coefficients of [II] 
and |III| should be independent of population, it follows 
that sated individuals are in no way influenced by the 
population of their own species or of prey, and to some 
degree can be said to be isolated from the community for 
the period of time they remain sated. It is reasonable 
to assume that for real animals subject to this condition, 
this period of time would be best spent in a den or other 
refuge offering relative safety from predation. Regard- 
less of the means by which such safety might be effected, 
which is more properly the concern of an agent-based 
model using spatial and individual-rationality considera- 
tions, we consider it proper to examine whether the vigi- 
lance of sated individuals has significant consequences for 
the model dynamics. We therefore denote as variant a of 
the model the case in which sated individuals are equally 
susceptible to predation as foraging individuals, and as 
variant b the case in which sated individuals cannot be 
predated at all. If it can be demonstrated that no signif- 
icant differences exist between these two cases, it is to be 
expected that reasonable intermediate cases can also be 
unproblematically related to the parent model. In so far 
as it makes no difference to the dynamics, we prefer to 
adopt variant b as the default case, since this suggests a 
more optimal individual behaviour. 
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In the parent model, foraging strategy is taken to adapt 
to the current population of all species on the shortest 
possible time-scale, so that at all times all individuals of 
all species follow the evolutionarily stable strategy, which 
is computationally arrived at by considering the effects 
on |3|) of prey availability and competitor abundance. 
A more plausible model of strategy selection would in- 
corporate at least the effects of limited knowledge, since 
individual predators must make decisions based on obser- 
vation or communicated information rather than directly 
based on global population data. Such decision making 
is also likely to be based on spatial distribution, and we 
therefore defer both informational and spatial consider- 
ations to the agent-based model we foresee as the next 
development of Webworld toward an ecosystem model 
grounded in individual-level considerations. For the pur- 
poses of this paper we consider a method of strategy se- 
lection such that the eventual, detailed behaviour model 
is likely to lie within the range investigated here. Vari- 
ant I denotes the case in which each time the population 
of any species changes, the strategy of each species is 
updated, and all individuals follow the ESS. For variant 
II we attempt to recover the ESS through population 
dynamics. In this case we assign to each individual a 
prey species which is their only prey for life, and are able 
to approximately track the ESS by causing offspring to 
adopt the same foraging strategy as their parent. Strate- 
gies which are under-represented in the population of a 
species have a higher reproduction rate than strategies 
which are over-represented, leading to their replenish- 
ment and convergence toward the ESS. A complication 
is that the population of each strategy is subject to de- 
mographic stochasticity, and hence at risk of extinction. 
With probability \x we assign each newborn individual 
to a strategy selected at random using the contemporary 
ESS as a probability distribution. This allows extinct 
strategies to be repopulated without losing large num- 
bers of individuals to unviable strategies. We investigate 
the results for /i in the range 0.01 < fi < 1, and use 
(i = 0.1 as the standard value. 



F. Summary 

The implementation of the stochastic Webworld model 
requires that each species, population Ni, be divided into 
sub-populations according to the number of prey it pos- 
sesses. The number of individuals actively foraging is 
given by Ny while the sated individuals are divided into 
groups according to the last prey item consumed, with 
populations Ni*j. Sated individuals return to the forag- 
ing population at a rate dependent on the prey species 
consumed, returning more rapidly for prey species corre- 
sponding to a large score, iSy. In variant I foraging effort 
is divided in the same manner as in the parent model, ac- 
cording to the evolutionarily stable strategy (ESS) given 
by (j20|) . In the more stochastic variant II each individual 
pursues a single prey species, and the relative reproduc- 



tive success of different foraging strategies contributes 
to the identification of the ESS by natural selection. For 
small forager populations it is necessary to avoid the per- 
manent extinction of strategies by allowing mutant off- 
spring to follow a strategy different from that of their 
parent, which occurs with probability /x. 

The second variation in the stochastic model that we 
consider is the ability to feed on sated individuals. It is 
implicit in the deterministic population dynamics that all 
individuals of a species are equally vulnerable to preda- 
tion, but as shown in Section III C| sated individuals are 
not influenced by any other individuals of their own or 
prey species. Following such work on optimal strategies 
as Houston et al. [HI, it mav be reasonable to assume 
that individuals not engaged in foraging are able cither 
to seek refuge where they are not vulnerable to predators 
or to adopt some vigilance strategy that removes the risk 
of predation at the expense of being unable to simulta- 
neously forage. This being the case, and noting that at 
the steady state the fraction of sated individuals is 

= , (34) 

N i3 bNj +J2k a ik fkj Skj N k ' 

which is typically of the order 10 -2 for N pley / N pTe d a t r ~ 
10 and S ~ 10, it is not strongly disruptive to the model 
to assume that only hungry individuals are susceptible 
to predation. Because we believe it to have the most 
plausibility, and most potential for developing a more 
detailed model, we use variant 116 as the default model 
against which others are compared. 

III. RESULTS BY PHENOMENOLOGY 

In order to compare the newly developed stochastic 
model with the parent model as closely as possible, we 
contrast their behaviour in three sets of circumstances. 
Each of these utilises the same 'mainland' community 
of species, evolved in the parent model and containing 
103 species, in order that incidental differences are min- 
imised. The simplest comparison between the parent and 
stochastic models is to directly compare the population 
dynamics, for which purpose it is preferable to make use 
of a food web of fewer species and lesser complexity than 
the mainland. Powell & McKane [22j introduced the 
method of using an evolved Webworld ecosystem as a 
species pool from which an 'island' community could be 
assembled by the addition of many species under the ac- 
tion of the population dynamics. Two small webs which 
resulted from this process, which we refer to as Islands A 
and B, are drawn in Table U and were used to study the 
population dynamics in the absence of the introduction 
of any further species. These results were repeated for 
several other webs of differing size as a test of generality. 
The labels associated with species indicate their trophic 
level, so species PI is a 'plant', HI a 'herbivore', and CI 
a 'carnivore'. Species E is the environment, population 
i?, which provides a constant food source for all plants. 
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Island A consists of two 
herbivore each feeding 
predominantly on pairs 
of plants. However, the 
most populous herbivore 
focuses some feeding 
attention on plant PI 
and also directly on 
herbivore HI 




Island B is already a 
highly complex web with 
interesting adaption to 
increasing stochasticity. 
The omnivory of the 
carnivore and herbivore 
H2 are not seen in the 
parent model. 




FIG. 1: The ensemble mean of time series generated with 
the stochastic population dynamics for Island A, using variant 
116 with /i = 0.1 and scaling £1 — 300. Population sizes are 
divided by Q for comparability. Solid lines mark species on 
the lower trophic level ('plants'); dashed lines mark species 
on the upper trophic level ('herbivores'). 



TABLE I: The food webs used in the comparison of time- 
series between the stochastic and parent population dynam- 
ics. In all webs, solid lines remain over all scalings. Long, 
medium and short dashed lines are those which are intermit- 
tent at f2 = 10, Q = 30 and Q — 100, respectively. The green 
arrows show intermittent links which exist only in the stochas- 
tic dynamics and are not represented in the parent model at 
all. The species designations are given in ascending order of 
population in the parent model. 



The second comparison of the two models is to create 
communities directly in the stochastic model. Because is- 
land creation restricts the set of species to those present 
in the mainland community, the results are more directly 
comparable than are those for independently evolved 
communities, and the process of immigration produces 
mature communities much more rapidly than is achieved 
by evolution from a single ancestral species. As such, the 
second test of the similarity of the two models was per- 
formed by statistical comparison of islands constructed 
under the two types of population dynamics. A third 
test, referred to as the 'reduction' process, is to take a 
large food web and record its behaviour as the resources 
supplied to basal species, R, are reduced. Using the 
mainland community as the initial food web, reduction 
of R tends to reduce the population of each species by 
the same factor until one or more species go extinct. Be- 
cause the food web structure changes at this point, the 
subsequent behaviour of the food web is not trivially pre- 
dictable. Under the population dynamics of the parent 
model the reduction process follows a deterministic tra- 
jectory, which can be compared to an example trajectory 
derived from the stochastic model. 

A simple feature of Webworld, which arises from the 
use of a ratio-dependent functional response, is that the 
dynamics of the parent model are completely unaltered if 
the basal resources, R, the population at which species go 



extinct and the population of each species are all scaled 
by the same factor. Equations @ and © show that, 
once scaling is taken into account, the population dy- 
namics of the system is unchanged. In the parent model 
it is therefore reasonable to set the scale by assigning the 
extinction threshold to a population of unity, but this 
does not necessarily correspond to one individual of the 
stochastic model. When comparing simple population 
dynamics with the parent model, we therefore adopt a 
scaling factor, Q, by which populations from the parent 
model are multiplied to obtain starting populations for 
the stochastic model. The results presented in Figs. [T] 
and [TO] are divided by the same scaling factor to 
assist comparison. The largest value of scaling adopted 
is fl = 300, which is a limitation of computational re- 
sources. The smallest value, = 10, is a consequence of 
the large fraction of simulation runs in which at least one 
species goes extinct within the first one hundred units of 
time, and suggests that the minimum viable population 
on this time scale is of the order ten. When considering 
island or reduction results, it is no longer appropriate to 
refer to a scaling, since species are removed from the com- 
munity by extinction by fluctuations about an unknown 
typical minimum population. Instead, Figs. 151 and [TT1 use 
the number of species present as a proxy for food web 
size when comparing other food web attributes between 
the parent and stochastic Webworld models. 

Having introduced variant II, in which the foraging 
strategy adopted by a species is not the evolutionarily 
stable strategy (ESS) but rather a collective strategy 
based on each individual feeding on a single prey species, 
it is important to assess the effect of choosing different 
values of the mutation probability, /i. The demographic 
extinction of certain strategies is inevitable if inheritance 
is important, and even when fi = 1.0 individuals retain 
a single strategy for their entire life-span, restricting the 
ability of the strategy of the species as a whole to adapt 
to the contemporary ESS. We therefore consider the ef- 
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FIG. 2: Selected species from Island A showing the effect of 
reducing the scaling from Q, = 300 (solid) through = 100 
(dashed) and ft = 30 (dot-dashed) to fi = 10 (dotted). Note 
that the y-axes are independently scaled and also overlap with 
each other. 

fects of five values of fi, spaced logarithmically in the 
range 0.01 < \i < 1.0, to which we add variant 16 as a 
comparable variant in which the ESS is followed at all 
times. 

The following three sections concentrate on three phe- 
nomena by which the parent and stochastic models can 
be compared. In Section IlII Al we examine how the par- 
ent model differs from the ensemble average of results 
from the stochastic model. In Section fill Bl we examine 
the nature of the fluctuations about the ensemble average 
seen in any given realisation of the stochastic model, and 
use these insights to understand some attributes of the 
webs persistent under the stochastic dynamics. In Sec- 
tion llll Cl we consider in detail the effects of using variant 
II &, in which links between species are formed and broken 
through the population demographics. 



A. Recovery of the parent model 

The ensemble-average time series from the stochastic 
dynamics in Island A at f2 = 300 and /i = 0.1 is repro- 
duced in Fig. [TJ Apart from the least populous plant, 
the ensemble average population of each species reaches 
a fixed point by the mid-point of the time-series which 
differs by less than 3% from that of the parent model, 
which is the initial condition. The discrepancy relates to 
the extent to which variant 116 is able to sustain those 
weak links which have been shown by McCann [25| to 
have an important influence on food web structure. No 
extinctions were observed for this value of fi, and ex- 
cellent agreement with the parent model is achieved for 
modest population sizes - each species present has a pop- 
ulation of less than ten thousand individuals. 

For smaller values of f2, the influence of demographic 
stochasticity is more pronounced. Not only do weak 
links become more susceptible to disruption, but species 
as a whole are more vulnerable to extinction. In or- 
der to understand a homogeneous data set we average 
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FIG. 3: The mean number of species for an ensemble of one 
hundred islands, against resources. Solid line; parent model. 
Dashed line; variant lib with /j, — 0.1. Cross; variant Ila with 
y, = 0.1. Circle; variant 16. 



over only those simulation runs in which extinction did 
not occur. Ensemble averages of the population dynam- 
ics, corresponding to Fig. [IJ are shown in Fig. [2] for 
the three species which showed significant deviation from 
the steady-state population found in the parent model. 
The probable cause of the deviations from the parent 
model result is the change in the food web structure as 
weak links are removed, and is discussed in detail in Sec- 
tion fllTC] The qualitative effect is that as smaller values 
of f2 are considered, the parent model fails to reflect even 
the ensemble average of the population dynamics. The 
important influence of weak links on the food web struc- 
ture is intimately connected with the population dynam- 
ics, since a weak link implies an event of rare occurrence 
whether or not the species involved are of small popula- 
tion. The ensemble average for small can be recovered 
in the parent model by specific selection of which links 
need to be incorporated, but this process cannot be prop- 
erly understood without first assessing the results of the 
stochastic model. 

The corresponding result for the island model is shown 
in Fig. [3J in which the number of species in a mature is- 
land community is shown for both the parent model and 
variant 116 of the stochastic model. The first point of 
agreement with the population dynamics results is that 
there exists an approximate factor of ten between the re- 
sources necessary to produce comparable results in the 
two models; if the curve for the stochastic dynamics is 
shifted a decade to the left, it approximately overlays 
the parent model results. In particular it is notable 
that complex communities are still readily formed by 
immigration under stochastic population dynamics, the 
largest communities constructed being limited by com- 
putational resources rather than intrinsic properties of 
the model. Noting that these results cor resp ond to the 
Type IV species-area curve of Scheiner [26j, which re- 
lates to the number of species in distinct islands and is 
therefore not necessarily monotonic, there is a good over- 
all agreement about the increase in species number with 
resource availability, which we equate to the physical ex- 
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FIG. 4: The mean number of species for an ensemble of one 
hundred islands against mutation probability fx, for R = 10 4 . 
The dashed line connects points belonging to variant 116; the 
final, disconnected point corresponds to variant 16 and hence 
is not part of the same interpolation. 



tent of a real ecosystem. The most significant deviation 
corresponds to communities at the smallest values of re- 
sources, where there is a step-like transition in the case 
of the parent model from systems with one or two species 
to those with several. This phenomenon was discussed in 
Powell & McKane [22 |. and is caused by the weakening 
of niche-sharing restrictions on plants when one or more 
herbivores enter the system. In the case of the stochastic 
population dynamics there is no noticeable step, and ex- 
amination of numerous example food webs indicates that 
even when herbivores become viable, it remains difficult 
to sustain multiple plant species. The underlying cause 
of the continued exclusion under stochastic dynamics is 
explored in Section lillBI The circle in Fig. [3] marks the 
number of species present when variant 16 is used rather 
than variant 116, and there is consequently no selection 
against weak links. The result is that fewer species are 
able to co-exist, a consequence of the improved ability of 
herbivores to exploit relatively rare plant species. Fig. 0] 
shows how the number of species simultaneously sus- 
tained for this value of resources, R = 10 4 , varies with the 
ability of species to adopt new strategies though chang- 
ing mutation, fi. The point corresponding to variant 16 
does not properly belong on the same scale as different 
values of /i, but conceptually corresponds to a situation 
of very large mutation rate and is therefore plotted to 
the right of the figure. It is notable that there exists a 
peak in the number of species, where small values of /x 
do not allow sufficiently diverse strategies to maximise 
the species diversity, but large values of /i also reduce 
diversity through the ease of resource exploitation, indi- 
cating a positive relationship between specialisation and 
diversity for this region of parameter space. 

A further check on the ability of the stochastic pop- 
ulation dynamics to support complex food webs is sup- 
plied by the reduction scenario. Fig. [5] shows the path 
that a single stochastic realisation of the reduction of 
resources follows, along with the unique path followed 
by the parent model. Similarly to the island results, 
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FIG. 5: The number of species remaining as resources are 
reduced, against resources. Solid line; parent model with slow 
reduction. Dashed line; variant 116 with /i = 0.1, reducing 
resources by 8% in unit time, compared with d — 1. 

there is approximately a factor of ten in resources be- 
tween the two curves, although this is partially obscured 
by the rapid rate at which it was necessary to constrict 
resources for the stochastic model. Whereas the popula- 
tion dynamics of the parent model can be run such that 
resources are reduced at an effectively infinitesimal rate, 
computational requirements and practical considerations 
require the stochastic model to be run at a definite re- 
duction rate. Computing the population dynamics of any 
community under stochastic population dynamics for an 
infinite length of time without speciation will result in 
the extinction of every species through fluctuations, so 
there does not exist an asymptotic curve for the stochas- 
tic model. Furthermore, the slowest rate of restriction 
for which results could be simulated corresponds to an 
8% reduction in resources per unit time period, which 
is rapid given that the expected lifetime of an individ- 
ual is also unit time. Species are therefore expected to 
persist to smaller values of R under these circumstances 
than would be true for a slower reduction rate, and the 
curve is somewhat smoothed by the stochastic delay in- 
troduced between the time a species becomes unviable 
and the moment at which it becomes extinct. Allowing 
for this effect, there is good agreement between the two 
curves, and the stochastic model is shown to be successful 
at allowing large, complex communities to persist. 

B. Fluctuations and intra-trophic neutrality 

One way in which the stochastic model differs from 
the parent model is that fluctuations in population sizes 
exist even when the ensemble average populations are 
constant. Extinction occurs when the fluctuating pop- 
ulation of a species drops instantaneously to zero, and 
hence the relative magnitude of fluctuations is an impor- 
tant contributing factor to food web structure. In this 
section we examine the variance of those fluctuations, and 
the extent to which the fluctuations in different species 
are not independent. Fig. [5] shows that, even for scaling 
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FIG. 6: The time series of one realisation of the stochastic 
population dynamics for Island A, using variant life with fj, = 
0.1 and scaling SI = 300. Blue lines; plant species. Black 
lines; herbivores. 



£1 = 300, these fluctuations are of significant amplitude, 
and it is notable that despite the similarity in popula- 
tion size, plant PI exhibits significantly larger fluctua- 
tions than herbivore H2. Fig. [7] shows the variance of 
each of the species averaged across an ensemble for each 
of the scaling factors, f2, examined. In every case the 
variance for plant species is larger than that for the her- 
bivores, the more populous of which has a variance nearly 
equal to its mean population for all values of used. For 
the larger values of the variance is approximately pro- 
portional to the population for all species, and hence the 
lines in Fig. [7] are approximately horizontal. This indi- 
cates that as the scaling factor increases, the standard 
deviation of population size increases only as the square 
root of the population. The ability to neglect fluctuations 
for very large population sizes is an assumption implicit 
in the parent model, and these results demonstrate that 
the fluctuations do indeed become negligible. 

If time series subject to independent noise are added, 
the variance of the summed time series is equal to the 
sum of the variance of each time series separately. The 
relative variance would therefore lie between the values 
calculated for each time series, yet Fig. [7J clearly shows 
that this is not the case for the stochastic Webworld 
results. Indeed, the variance for the sum of the plant 
populations is nearly equal to the total population, as 
is also the case for the sum of herbivore populations. 
This indicates that the fluctuations in the plant popu- 
lations must be anti-correlated. In Fig. [8] a time series 
for plant P4 is shown along with that for the sum of the 
other three plants. It can be seen that significant anti- 
correlation does exist, which is especially obvious for the 
large fluctuations around time 55 and 70. To quantify 
this relationship we calculate the correlation coefficient, 
r 2 , between time series [27j . It was found that an in- 
crease in r 2 could be found by allowing a offset between 
the time series of different species, and Table |TT] records 
the maximum value of r 2 and the corresponding time off- 
set, along with the gradient of the best-fitting regression 
line. The results show that for Island A two pairs of 
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FIG. 7: The ratio of the variance of population time series to 
the mean population at various values of scaling, fi. Dashed 
lines; plant species. Dotted lines; herbivore species. Solid 
line; total plant population. Dash-dotted line; total herbivore 
population. 



plants exist with relatively strong correlation; species PI 
correlating strongly with P3, and P2 with P4. As shown 
in Table HI these pairs form trophic species in the full 
sense that they share both food source and predators, al- 
though PI is also weakly predated by H2. These trophic 
species are strongly correlated with each other, and are 
more strongly correlated their respective herbivores than 
is any individual plant species. The optimal time offset 
for any pair of species within the same trophic level is 
near zero, while the best correlation between the trophic 
plants and herbivores is found with the herbivore time 
series lagging by up to approximately unit time, suggest- 
ing that changes in the herbivore population can, in part, 
be attributed to fluctuations in the abundance of their 
food source. The gradient of the regression line between 
each pair of plant species is close to — 1. This is a conse- 
quence of the strong competition between species which 
share a common food source. Although the total popula- 
tion across all plant species is relatively well-constrained, 
any one plant species can increase in population at the 
expense of it competitors. Meanwhile there is a positive 
correlation between the population of each herbivore and 
that of the trophic species upon which it predominantly 
feeds. With the opposite gradient, there is a consequent 
correlation between herbivore population and the com- 
petitors of their prey. 

A similar analysis was carried out with Island B, of 
which the results for the sub- web comprising PI, P2, HI 
and H3 are particularly notable and are appended to Ta- 
ble ITT] Two major differences between this sub- web and 
Island A are the presence of stronger cross-links and the 
presence of the carnivore, making the combination of HI 
and H3 a trophic species. As well as the correlation be- 
tween plants, there is a significant correlation between 
herbivores, also with near-zero offset. However, although 
the optimal offset between herbivore and plant time se- 
ries is of similar magnitude to those seen Island A, in 
Island B it is the herbivore populations which lead. The 
gradient in Table |TT] is equal to the number of individ- 
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TABLE II: Selected data from the analysis of correlations in 
both Islands A and B. The j/-time series is offset by At with 
respect to the :r-series to maximise the correlation coefficient, 
r 2 . The gradient shown is that of the best-fitting regression 
line at At. Note that the typical population of species in 
Island B is greater than in Island A, contributing to larger 
values for r 2 . 



uals of the plant species necessary to support a single 
extra herbivore. This number is greater in Island B due 
to the additional losses to predation by the carnivore, 
whose inclusion may be responsible for altering the rela- 
tive strength of top-down and bottom- up effects, reflected 
in the change from lagging to leading herbivore fluctua- 
tions. 

In the parent model, the co-existence of multiple 
species within a trophic species is allowed because inter- 
specific competition is less that intra-specific competi- 
tion, encoded in the competition factor, a. This can be 
interpreted as reflecting differences in the feeding habits 
of competitor species which are not explicitly modelled. 
Under stochastic population dynamics an additional ef- 
fect is observed in which demographic fluctuations in 
the populations of the niche-sharing species are larger 
than those of species which do not share a niche. By 
extension we expect that an increasing overlap of com- 
mon predators and prey between two species corresponds 
to an increasing amplitude of fluctuations in population, 
while the total population of the niche (trophic species) 
remains relatively steady. As inter-specific competition 
increases, and a — » 1, both the parent and stochastic 
models would predict the extinction of the less suitable 
taxonomic species within trophic species. 

A consequence of the large fluctuations in the popula- 
tion of trophically similar species in the stochastic model 
is a stronger exclusion of niche-sharing species. In turn, 
fewer plants are able to co-exist in island communities. 
This is reflected in the reduced number of species present 
in stochastic islands relative to parent-model islands, as 
can be seen in Fig. [3] A further effect, shown by Fig. [51 is 
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FIG. 8: An example of anti-correlations between populations 
of basal species, taken from Island A using variant 116 with 
jit = 0.1, scaling Q — 300. The lower line marks the population 
of the most populous plant, P4. The upper line marks the 
total population of the other three plants. 

that the average number of trophic levels is higher in the 
case of the stochastic model than for the parent model 
when comparing webs with an equal number of species 
and no carnivores (approximately, S < 10). Since webs 
containing multiple plants and no herbivores are less sta- 
ble under stochastic dynamics, it follows that for a given 
number of species there will be more webs with at least 
one herbivore. 



C. Strategy extinction and web restructuring 

Having introduced demographic effects into foraging 
strategy through variant II of the stochastic model, in 
which each individual feeds on a single prey species for 
its entire life-span, species become susceptible to 'forget- 
ting' feeding links if the typical number of individuals 
following that strategy is small. As the scaling factor, 
f2, is reduced, species become increasingly susceptible to 
this effect, and departures from the food web structure of 
the parent model are expected to increase in magnitude. 
The diagrams of the Island A and B food webs shown 
in Table [J illustrate how weaker feeding links are succes- 
sively lost as Q is decreased from 300, when the webs 
are essentially equivalent to those in the parent model, 
through ft = 100, 30, and 10. Once a feeding link is 
removed from the model, the prey species typically in- 
creases in population and the predator decreases, since 
it is receiving income from fewer prey individuals. The 
context within the food web structure complicates this 
picture, since the removal of competition strengthens the 
ability of other predators to exploit a species. For exam- 
ple, if the link between PI and H2 in Island A is lost, HI 
does not suffer inter-specific competition when feeding on 
PI, and is therefore able to increase in population. This 
reduces the beneficial effect to PI of this link being lost. 
A second example of the indirect consequences of the loss 
of weak links is provided by the cessation of feeding by 
H2 on HI. When this occurs HI increases in population, 
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FIG. 9: The mean number of trophic levels against the mean 
number of species present, for an ensemble of one hundred 
islands at each point. Solid line; parent model. Dashed line; 
variant 116 with [i — 0.1. Cross; variant Ila with [i — 0.1. 
Circle; variant 16. 



FIG. 10: Parametric plot of the population of plant PI and 
the two herbivores of Island A, using variant 116 with scaling 
Q = 50. From right to left the points correspond to /i = 0.01, 
0.03, 0.1, 0.3, and 1.0. The disconnected points to the left 
correspond to variant 16. 



increasing the predation pressure on plants PI and P3, 
whose populations decrease. The propagation of complex 
indirect effects through the food web can be understood 
in basic terms by investigation of the proximal effects 
of link removal, but the net effects on the community 
cannot be predicted any more simply than by directly 
simulating the population dynamics, due to the complex 
nature of competition and predation relationships. 

In addition to the loss of weak links, some links appear 
in the stochastic model which would be predicted from 
the parent model to correspond to less than a single in- 
dividual. For small values of f2 the food web structure 
has been changed by the deletion of other links, with the 
result that the populations of species differ from their val- 
ues in the parent model. Under these circumstances links 
may become viable for one of two reasons; either the prey 
species has increased in abundance to the point at which 
it forms an exploitable resource, or competition for that 
prey has decreased. In general both of these effects oc- 
cur, and several such links are marked on the diagram of 
Island B in Table HI All these links are intermittent - the 
predator species does not permanently have individuals 
which follow the strategy, but that sub-population goes 
extinct and is later re-established through the birth of a 
mutant individual. In the case of extremely simple food 
webs it can be seen that there are population changes in 
the prey species depending on whether or not it is be- 
ing fed upon, but in even moderately complex food webs 
such effects are within the demographic noise and hence 
cannot be clearly demonstrated. 

Fig. ITTJI demonstrates how the population of plant PI 
increases with decreasing fi, which corresponds to an in- 
crease in the expected time before an individual of species 
H2 mutates to re-populate the extinct foraging strategy. 
Correspondingly, herbivore H2 is seen to decrease in pop- 
ulation with decreasing /i, since this source of nutrients 
forms an increasingly small part of its diet. Herbivore 
HI, whose feeding link to PI does not significantly de- 
crease over this range of fi, experiences a population in- 



crease correlated with the increasing abundance of its 
prey. Because HI does not feed exclusively on this plant, 
its population increase is smaller as a fraction of its total 
population than the increase in PI, which becomes an 
increasing fraction of its diet with increasing abundance. 
This figure demonstrates specific causes of both positive 
and negative correlations between populations, illustrat- 
ing the complex outcomes of the food web dynamics. 

The large-.!? part of Fig. [9] shows that the number 
of carnivores is less in the stochastic than the parent 
model when considering food webs with equal numbers of 
species, but that carnivores are not completely excluded 
from the food webs. The reduction in their prevalence is a 
consequence of the removal of weaker feeding links, since 
carnivores in the Webworld model tend to have more di- 
verse diets than do lower trophic levels. As an example 
of this, the single carnivore in Island B feeds on each of 
the four herbivores, but no herbivore in the parent model 
(i.e. excluding shaded links in Table fl} feeds on more than 
two species. Because the population required to sustain 
links to four prey species is greater than that necessary to 
feed on a single prey, the entry of carnivores into the food 
web is delayed relative to the entry of the lower trophic 
levels, and in particular it is possible to support a larger 
number of plant and herbivore species before carnivorous 
strategies become viable. As a consequence, the curve for 
the stochastic model lies below that for the parent model 
when considering trophic levels greater than two. The 
absence of these feeding links also increases the fraction 
of 'top' (unpredated) species, as shown in Fig. [TTJ For a 
small total number of species, S, the curves for the parent 
and stochastic models are very similar as only plants and 
herbivores exist, and despite differences in the proportion 
of species found in those two trophic levels as discussed 
in Section fill Bl The circle marks the point correspond- 
ing to variant I&, in which the presence of feeding links 
is not affected by demographic effects, and which as a 
consequence is within the region occupied by the parent 
model. The inset to Fig. QT] shows the monotonic de- 
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FIG. 11: The mean fraction of unpredated species against the 
mean number of species present, for an ensemble of one hun- 
dred islands at each point. Solid line; parent model. Dashed 
line; variant 116 with fi = 0.1. Cross; variant Ila with fi = 0.1. 
Circle; variant 16. Inset: The mean fraction of unpredated 
species against mutation probability (i, for R — 10 4 . 

crease in the fraction of top species as \x is increased, and 
the continuation of this trend into variant 16. 



D. Refuge for sated individuals 

The distinction between variants a and 6 of the model 
corresponds to the ability of sated individuals, who can- 
not gain advantage from any interaction, to avoid inter- 
action with predators. In variant 6 sated individuals are 
completely safe from predation, whereas in variant a they 
are as vulnerable as foraging individuals are. The plau- 
sible behaviour of a more detailed model is likely to lie 
between these two extremes. Table lllll records the ensem- 
ble average population for each of the six species of Island 
A for variants Ila and lib, and variant 16 for compari- 
son. The standard deviations, also shown, indicate the 
size of the population fluctuations. It is apparent that for 
each species the difference between variants Ila and 116 
is much smaller than the corresponding fluctuations, but 
the differences between variants I and II are significant. 
These findings are reflected in the position of markers 
in Figs. [31 O and [TT] In each case variant Ila produces 
results essentially equivalent to variant 116, while variant 
16 shows significant differences. In the latter two figures, 
where the number of species in an island food web is used 
to relate the size of parent-model and stochastic-model 
webs, variant 16 produces results similar to the parent 
model. If refuge has important consequences for food 
web structure then its effects arc likely to be concealed 
within the functional response used by Webworld. 

IV. CONCLUSIONS 

The primary aim of this paper was to develop the Web- 
world model by introducing an individual-based popu- 
lation dynamics. The Webworld model has been sue- 





Variant Ila 


Variant 16 


Variant 116 


Species 


M 


a 




a 


/i a 


PI 


12.78 


4.8 


5.9 


2.7 


13.10 4.5 


P2 


15.49 


3.3 


16.8 


2.9 


15.64 3.2 


P3 


18.44 


3.6 


21.6 


3.0 


18.51 3.6 


P4 


24.07 


3.1 


25.6 


3.0 


24.24 3.1 


HI 


2.31 


0.86 


1.28 


0.41 


2.326 0.85 


H2 


4.22 


0.82 


4.60 


0.72 


4.164 0.81 



TABLE III: A comparison of stationary behaviour between 
variant Ila, in which all individuals are vulnerable to preda- 
tion, and variant 116, in which only hungry individuals are 
vulnerable, is shown for Island A. The scaling used here is 
f2 = 10. fi; mean, a; standard deviation. 



cessful in constructing food webs with several realistic 
properties, based on simple population dynamics and 
foraging strategy optimisation. By developing a model 
whose basic assumptions correspond to a finer level of 
detail, the interactions between individuals rather than 
species-level descriptions, it becomes possible to explore 
the consequences of a range of behaviours enacted at the 
individual-level. For instance, at the level of the par- 
ent model it is possible to test robustness to different 
forms of the functional response, whose normal form is 
given by ([3]), but use of an individual-based model al- 
lows the more direct comparison of different forms of 
the reaction rate equations which are assumed to un- 
derlie species-level interactions. In addition, it has been 
shown that the stochastic model closely approximates the 
parent Webworld model when population sizes are large, 
but that when considering constructed food webs this is 
never an acceptable assumption, since there always exist 
some species which would not be viable for a significantly 
smaller ecosystem, and whose population demographics 
have a profound influence on the community as a whole. 

A secondary aim in developing an individual-based ver- 
sion of the Webworld model was to establish a point on 
the way to an agent-based model. At the level of the indi- 
vidual, it becomes possible to model ecologically impor- 
tant processes such as predator-evasion, foraging choice 
through individual access to observations about prey dis- 
tribution, and to include the effects of life-history de- 
cisions and mating strategy. Each of these can be ex- 
pected to have effects on the population dynamics of 
the species and hence on the structure of the ecosystem, 
but it is not practically possible to create a simulation 
model which takes individuals with detailed behaviour 
as a starting point and from that constructs a food web 
structure. Instead we have adopted the strategy of tak- 
ing a model known to produce realistic food webs, within 
which framework we can examine at various scales the in- 
fluences that differences in the details of one level of ab- 
straction have on the outcome of a broader picture. If a 
bridge can be built linking choice at the level of the indi- 
vidual at one end to ecosystem behaviour at the other, it 
will be possible to explore the consequences of the former 
on the latter despite the practical impossibility of com- 
puting the whole of an ecosystem at a sufficient level of 
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detail to resolve detailed variation between individuals. 

Some of the important results to come out of the pre- 
liminary work on the stochastic model performed in this 
paper relate to the effects of demographic stochasticity 
in different food web positions. In particular it has been 
shown that interactions between species significantly al- 
ter the magnitude of demographic fluctuations, and that 
taxonomic species forming part of a larger trophic species 
are more susceptible to these fluctuations than are species 
who share neither predators nor prey. The stochastic 
model broadly reproduces the results of the parent Web- 
world model in terms both of population dynamics and 
the construction of communities through immigration, 
although it is necessarily the case that extinction in the 
stochastic model occurs through fluctuations occasion- 
ally reducing a population to zero, rather than the mean 
population reaching a certain threshold as is the case in 



the parent model. Webs constructed in the parent model 
are robust to the imposition of the stochastic population 
dynamics, and a future line of enquiry would be to estab- 
lish whether this is due to the degree to which the two 
sets of population dynamics are equivalent, or whether 
the method of constructing food webs through an evolu- 
tionary process of species addition selects for robustness 
to changes in relative species abundance, and hence pro- 
vides stability against the demographic fluctuations. 



Acknowledgements 

We wish to thank EPSRC (UK) for funding; CRP un- 
der grant number GR/T11784 and RPB under a post- 
graduate grant. 



[1] Bulmer, M., 1994. Theoretical Evolutionary Ecology 
(Sinauer Associates Inc) 

[2] Blackburn, N., Azam, F., Hagstrom, A., 1997. Spatially 
explicit simulations of a microbial food web. Limnology 
and Oceanography 42 613 

[3] Vatland, S., Budy, P., Thiede, G. P., 2008. A bioenerget- 
ics approach to modeling striped bass and threadfin shad 
predator-prey dynamics in Lake Powell, Utah-Arizona. 
Transactions of the American Fisheries Society 137 262 

[4] Williams, R. J., Martinez, N. D., 2000. Simple rules yield 
complex food webs. Nature 404 180 

[5] Solow, A. R., Beet, A. R., 1998. On lumping species in 
food webs. Ecology 79 2013 

[6] Drossel, B., Higgs, P. G., McKane, A. J., 2001. The in- 
fluence of predator-prey population dynamics on the long 
term evolution of food web structure. J. Th. Biol. 208 91 

[7] Drossel, B., McKane, A. J., Quince, C, 2004. The im- 
pact of nonlinear functional responses on the long-term 
evolution of food web structure. J. Theor. Biol. 229 539 

[8] Quince, C, Higgs, P. G., McKane, A. J., 2005. Deleting 
species from model food webs. Oikos 110 283 

[9] Quince, G, Higgs, P. G., McKane, A. J., 2005. Topo- 
logical structure and interaction strengths in model food 
webs. Ecol. Model. 187 389 
[10] Vik, J. O., Brinch, C. N., Boutin, S., Stenseth, N. G, 
2008. Interlinking hare and lynx dynamics using a cen- 
tury's worth of annual data. Population Ecology 50 267 
[11] Saether, B. E., Engen, S., Islam, A., McCleery, R., Per- 
rins, G, 1998. Environmental Stochasticity and Extinc- 
tion Risk in a Population of a Small Songbird, the Great 
Tit. Am. Nat. 151 441 
[12] Lande, R., Engen, S., Sasther, B. E., 2003. Stochastic 
Population Dynamics in Ecology & Conservation (OUP) 
[13] Kfivan, V., Eisner, J., 2006. The effect of the Rolling type 
II functional response on apparent competition. Theoret- 
ical Population Biology 70 421 
[14] McKane, A. J., Newman, T. J., 2005. Predator-prey cy- 
cles from resonant amplification of demographic stochas- 



ticity. Phys. Rev. Lett. 94 218102 

[15] Magurran, A. E., Henderson, P. A., 2003. Explaining 
the excess of rare species in natural species abundance 
distributions. Nature 422 714 

[16] Powell, G R., McKane, A. J. Predicting the species abun- 
dance distribution using a model food web. 
larXiv:0805.0084t T 

[17] Lomolino, M. V., 2000. A call for a new paradigm of 
island biogeography. Global Ecology & Biogeography 9 1 

[18] Melbourne, B. A., Hastings, A., 2008. Extinction risk 
depends strongly on factors contributing to stochasticity. 
Nature 454 100 

[19] Gillespie, D. T., 1977. Exact stochastic simulation of cou- 
pled chemical-reactions. J. Physical Chemistry 81 2340 

[20] Lugo, C. A., McKane, A. J., 2008. The characteristics 
of species in an evolutionary food web model. J. Theor. 
Biol. 252 649 

[21] Lugo, G A., McKane, A. J., 2008. The robustness of the 

Webworid model to changes in its structure. Ecological 

Complexity 5 106 
[22] Powell, C. R., McKane, A. J. Comparison of food webs 

constructed by evolution and by immigration. 
larXiv:08 08.2922vl 
[23] Huisman, G., De Boer, R. J., 1997. A formal derivation 

of the "Beddington" functional response. J. Theor. Biol. 

185 389 

[24] Houston, A. I., McNamara, J. M., Hutchinson, J. M. G, 
1993. General results concerning the trade-off between 
gaining energy and avoiding predation. Phil. Trans. R. 
Soc. Lond. B 341 375 

[25] McCann, K. S., 2000. The diversity-stability debate. Na- 
ture 405 228 

[26] Scheiner, S. M., 2003. Six types of species-area curves. 
Global Ecology and Biogrography 12 441 

[27] Hogg, R. V., Craig, A., 1994. Introduction to Mathemat- 
ical Statistics (Prentice Hall), 5th edn. 



